Esaki spin diode based on composite graphene superlattices 
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Since the pioneering work by Esaki , quan- 
tum tunneling and negative differential resistance 
(NDR) have been the underlying principle of op- 
eration of various quantum devices . The prin- 
ciple implies that the electric current through a 
device increases when the chemical potential in 
one of the leads approaches an energy level of 
the device. However, if the energy level depends 
on the applied voltage the current can decrease 
drastically with a further increase of the voltage, 
resulting in the NDR. 

Such conductance anomaly can be observed 
in semiconductor heterostructures , molecular 
systems and at the atomic scale '. Here we ad- 
dress a novel device based on graphene, whose 
remarkable charge transport properties and long 
spin-coherence length makes it a promising ma- 
terial for spintronics • . We consider a compos- 
ite superlattice comprising a set of ferromagnetic 
insulator strips deposited on top of a graphene 
nanoribbon (GNR). We calculate current voltage 
characteristics of the device and predict spin- 
filtering effect and strong spin-dependent NDR. 

The proposed system is composed of a rectangular 
GNR of width W ~ 9.8 nm, connected to source and 
drain leads, and N = 5 rectangular strips of a ferro- 
magnetic insulator arranged periodically on top of the 
GNR (see the upper panel of Fig. 1). The width of the 
strips is a = 23.9 nm and the spacing between them is 
b = 55.8 nm. It is known that both the width and the 
edge type of a GNR strongly affect its electronic prop- 
erties. Both the experimental evidence 10 and ab-initio 
calculations show that the spectrum of a GNR with 
armchair edges has a band gap, which is inversely pro- 
portional to the width W and depends on the remainder 
(W/ao mod 3) (where ag = 0.246 nm is the width of 
the graphene lattice hexagon). Contrary to GNRs with 
zigzag edges, the dispersion relation of the armchair GNR 
is centered around k = 0; the latter is advantageous for 
tunneling structures, such as a superlattice, because res- 
onant levels are expected to be broader and less affected 
by disorder 12 . We therefore restrict ourselves to the arm- 
chair GNRs. 

EuO can be used as the ferromagnetic insulator; this 
material has been studied in conjunction with graphene 
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FIG. 1. Device design and potential profiles. Upper 
panel shows the schematics of the device: a GNR, connected 
to source (S) and drain (D) leads, with 5 perpendicular strips 
of a ferromagnetic insulator (green bars) placed on top of 
it. Potential profiles for spin-up (dotted red lines) and spin- 
down (dashed blue lines) electrons in the unbiased and biased 
devices are shown in the middle and lower panels respectively. 



both experimentally and theoretically 14 . The exchange 
interaction between magnetic ions in the strips and 
charge carriers in the GNR can be described as an ef- 
fective Zeeman splitting ±A 0X of the spin sublevels 1 . 
There is still no consensus on the magnitude of the ex- 
change splitting amplitude A ex in graphene. We use 
A ox = 5 meV, which lies in the range of values known 
from the literature (3 — 10 meV) 1 1-1,1 . We have checked 
that our results do not change qualitatively if we use a 
different value of A cx from the known range. 

Because the exchange interaction has the characteristic 
length scale of one atomic layer, the splitting is induced 
only in the regions of the GNR which are just below 
the ferromagnetic strips. Therefore, for the chosen sys- 
tem geometry, a spin-up (spin-down) electron propagat- 
ing along the sample will be subjected to a potential com- 
prising a periodic set of rectangular barriers (wells), as 
plotted in the middle panel of Fig. 1. In other words, the 
array of the ferromagnetic strips creates a spin-dependent 
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FIG. 2. Transmission and its polarization for the unbiased device, a, transmission coefficients T± as functions of 
energy for spin-up (red) and spin-down (blue) electrons. There is a very good agreement between the Dirac transfer-matrix 
theory (solid lines) and the tight-binding calculation (dotted lines). The origin of the energy for each curve is set to the 
lowest subband edge Eq calculated within the corresponding model. Horizontal bars in the lower part of the figure indicate 
transmission bands of the infinite superlattice. b, transmission polarization Pt as a function of energy for different numbers 
N of ferromagnetic strips. 



superlattice. We note that similar systems manifesting 
the NDR have been studied 1 ' , however the superlattice 
potential was supposed to be induced by electrostatic 
gates, so all characteristics were spin independent. 

Calculation methods and techniques. For not too nar- 
row GNRs the low energy excitations can be treated 
very efficiently within the Dirac approximation ■ . The 
boundary conditions at the GNR edge can be met by a 
superposition of two states from different valleys. The 
transverse momentum k± n (measured from the respec- 
tive Dirac point) is quantized"", and the longitudinal 
momentum fen obeys the dispersion relation E n (kn) = 
±HvF(k'j_ n + fcjj) 1 / 2 . In our case W — 40 a and the min- 
imum transverse momentum is k± = ir/[3(W + ag)} . 
The resulting gap is an important feature of the model 
because it eliminates the Klein tunneling effect which 
would always result in full transmission and current leak- 
age. 

For the chosen geometry, the potential depends only 
on the longitudinal coordinate x, and the transverse mo- 
mentum k± is conserved. Therefore it suffices to solve 
for the longitudinal wave function <j>\\(x). The solution 
of the Dirac equation within each interval of constant 
potential value Vj is the superposition of two counter- 
propagating sublattice pseudo-spinors. Both the longitu- 
dinal momentum kj and the spinor components depend 
on the energy E of the injected carrier and the local po- 
tential Vj. At each junction the wave vector kj changes 
while the wave function remains continuous. The free 
propagation across the intervals and the continuity at 
the interfaces can be cast in a transfer-matrix formalism, 
which ultimately yields the transmission and reflection of 
the (finite) system, as function of energy 21 . 

Transmission at zero bias. Because the superlattice 
potential depends on the carrier spin, the transmission 
coefficient T± is also spin-dependent. Hereafter, + (— ) 
signs and red (blue) colors in all figures correspond to 
spin-up (spin-down) electrons respectively. Figure 2 a 



shows the transmission probability through the unbiased 
sample calculated within the Dirac approximation (solid 
lines) and the full tight-binding model (dotted lines), 
which is more accurate" 1 . The figure demonstrates very 
good agreement between the two approaches. The Dirac 
model overestimates slightly the lowest subband edge po- 
sition Eq (by 0.3%), due to the k ■ p approximation in 
the derivation of the Dirac equation, which breaks down 
for very narrow GNRs~ . In the following, we use exclu- 
sively the Dirac transfer-matrix method since it requires 
less computational resources. Furthermore, the formal- 
ism can be applied to study the system as the number of 
strips N —> oo, and gain physical insight into the main 
features of the transmission 21 ' 22 . For an unbiased infinite 
superlattice, an eigenfunction has the Bloch phase factor 
exp(±iKl) which is the eigenvalue of the transfer matrix 
G across one superlattice cell of the length I = a + b. 
Thus, the dispersion relation for a GNR is implicitly de- 
fined by the equation cos(Kl) = Tr(G)/2. If | Tr(G)| > 2, 
there is no propagating solution and the corresponding 
energy E falls into the bandgap of the superlattice. 

The transmission spectrum, shown in Fig. 2 a, is typ- 
ical for a superlattice. Already for N = 5, the re- 
gions of high transmission coincide quite well with the 
bands of the infinite superlattice given by the horizontal 
bars at the bottom of the figure. The spectrum is spin- 
dependent, which manifests itself clearly in the transmis- 
sion polarization, defined as Pt = {T + — T_)/ (T + +T_), 
shown in Fig. 2 b. 

As the number of strips is increased, the transmission 
of modes with energies outside the transmission bands 
vanishes rapidly, thus leading to an enhanced polariza- 
tion. For TV > 3, the transmission polarization mani- 
fests rapid changes in almost the whole range of its pos- 
sible values (—1, f ) within narrow energy intervals. Such 
abrupt polarization switching can be expected only if the 
overlap between transmission bands corresponding to dif- 
ferent spins is small (see Fig. 2 a). The band overlap is 
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FIG. 3. Transmission and currents at finite bias. Panels a and b show that the transmission bands for both spins at 
finite bias Vsd are shifted, quenched and distorted compared to the unbiased case (Fig. 2 a). Panels c and d display the spin- 
polarized currents I± as functions of the bias Vsd, for T = 4K and different values of the chemical potential: fj, — Eo = OmeV 
(dashed), 0.5 meV (dotted), and 1.0 meV (dash-dotted) lines. Finally, panels e and f show the total current / and the current 
polarization P for the same parameters. 



determined by the relationship between different param- 
eters: the splitting amplitude A ex and the superlattice 
geometrical parameters a and b, which should be chosen 
carefully in order to observe a pronounced switching and 
filtering effect in a real world device 21 . Such a choice 
can be made, for example, by analyzing the transmission 
spectrum of the infinite lattice in the Dirac model. 

Spin-polarized current at finite bias. We finally use 
the obtained transmission coefficients T± to calculate 
the spin-polarized electric current within the Landauer- 
Biittiker scattering formalism 21 ' 2,5 . We address the cur- 
rent and its polarization at a finite bias Vsd between 
source and drain, whose chemical potentials, fis = M + 
eVsD and ^d = M; have the same offset [i from the band- 
gap center (the Dirac point). For simplicity, we assume 
that the source-drain voltage drops at the edges of the 
wells and barriers of the superlattices only, resulting in 
a piecewise potential profile as shown in the lower panel 
of Fig. 1. We have checked numerically that our results 
do not depend crucially on the details of the potential 
profile. The bias results in a distortion of the transmis- 
sion bands: the bands shift to lower energies, become 
quenched and finally disappear as the voltage increases, 
as shown in Fig. 3 a and b. 

The polarized currents I± as functions of Vsd are plot- 
ted in Fig. 3 c and d. The spin-dependent transmission 
bands and their distortion due to the bias lead to NDR 
regions at different values of the bias voltage for different 
spins. For the spin down, the NDR occurs at a lower bias 
and the negative slope of the I — V curve is particularly 
steep. This is due to the fact that the first transmission 
peak remains very sharp while it is being quenched (see 



panels a and b of Fig. 3). The lowest spin- up transmis- 
sion band disappears at a higher bias, when its profile is 
more washed out, resulting in the less pronounced NDR. 

It is worthwhile to address the total current / = 
I + + I- through the device, as well as its polarization 
P = (1+ — I-)/{I + + /_). These quantities are plotted 
in the panels e and f of Fig. 3. Panel e shows that the 
total current I also manifests the NDR for two different 
biases, corresponding to the NDR regions of /_ and I+. 
The current polarization shows an initial range with neg- 
ative values followed by a second region dominated by the 
spin- up current. As the bias increases further the polar- 
ization decays and finally vanishes. Note that the current 
is highly polarized for certain biases, which proves that 
the device can operate as a spin filter. On the other 
hand, because the characteristics I±(Vsd) are essentially 
different, if the source feeds partially polarized electrons, 
the total current through the device would depend on the 
degree of their polarization. The latter opens a possibil- 
ity to determine the polarization of a current by purely 
electrical measurements, which is a very promising appli- 
cation. 

We have considered an ideal device with perfect rect- 
angular GNR and strips, while different imperfections 
and perturbations can introduce disorder into the sys- 
tem and affect the electric current and its polarization 24 . 
There are different possible sources of disorder, for exam- 
ple, charged impurities in the substrate or defects of the 
device fabrication. The former results in an additional 
smooth electrostatic potential and can hardly deteriorate 
the transmission through the device to a large extent. 
However, the effect of the latter on the transport proper- 
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ties can be stronger. To estimate the possible impact of 
the fabrication imperfections on the predicted effects, we 
considered disordered superlattices varying randomly the 
strip widths and spacings by up to 20%. Our calculations 
demonstrate that the transmission bands are affected by 
the disorder to a comparable degree for both spin up 
and spin down electrons, which suggests that a moderate 
disorder would not seriously deteriorate transport and 
polarization properties of the device. The current mag- 
nitude remains almost the same, and the NDR turns out 
to be robust under the effects of disorder as well. 

In summary, we propose a novel graphene-based de- 
vice comprising a GNR and a regular array of ferromag- 
netic strips on top of it. The ferromagnets induce an 
exchange splitting of electronic states in the GNR and 
create a spin-dependent superlattice. We have shown 
that the electric current through the device can be highly 
polarized. Thus, the device can operate as a spin fil- 
ter. Alternatively, it can be used to obtain the po- 
larization degree of the source electrons (the polariza- 
tion of the "incoming" current) by purely electrical mea- 
surements. Moreover, the two polarized components of 
the current manifest non-monotonic dependencies on the 



source-drain voltage. In particular, for both spins, the 
current-voltage characteristics presents regions with neg- 
ative differential resistance for the bias in the range of a 
few millivolts. The device operates therefore as a low- 
voltage Esaki diode for spin-polarized currents. 

A prominent advantage of the usage of a superlattice 
induced by ferromagnets is that the exchange interaction 
is very short-ranged; its characteristic length scale is on 
the order of one monolayer. Unlike the long-range elec- 
trostatic gate potentials which can interfere with each 
other, setting a practical lower limit for the inter-device 
spacing, the exchange interaction induced potential pro- 
files are very abrupt. Therefore, heterostructures created 
by ferromagnets allow for very close packing of circuits 
and, consequently, considerably higher device densities. 

Finally, we note that in a spintronic device the degree 
of freedom that carries information is the polarization of 
the current rather than its magnitude. We have shown 
that the current polarization is also a non-monotonic 
function of the source-drain voltage in our proposed de- 
vice, which suggests that it can be used as an Esaki spin 
diode. This opens a possibility to design a whole new 
class of true spintronic circuits such as spin oscillators, 
amplifiers and triggers. 
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SUPPLEMENT A. MODEL HAMILTONIAN 



A simple tight-binding (TB) Hamiltonian of a single electron is widely used to model GNRs. For 
low energy excitations, i.e., energies close to the Dirac point, the interaction can be restricted to 
nearest neighbors. Then, the Hamiltonian can be written as 



Here |z) is the ket vector of the atomic orbital of the ith carbon atom, t = 2.8 eV is the hopping 
between neighboring atoms, the full set of which is denoted as The on-site energy is the 

superposition of the following two terms: the bias induced electrostatic potential ego(xi), where x, 
is the coordinate of the ith atom along the direction of the GNR, and the spin- dependent exchange- 
interaction shift A ex due to the ferromagnetic strips, with o = ±1 for spin-up and spin-down 
electrons. The exchange-interaction is induced only at atoms that are in direct contact with the 
ferromagnetic strips (the full set of them is labeled as C in the above equation). Both terms are 
sketched in the middle and lower panels of Fig. 1 of the main text. 

First, we calculate the wave function in the whole sample and the transmission coefficient. This is 
accomplished by assuming semi-infinite leads, whose modes are calculated using an effective transfer- 
matrix approach 25 . Then, both the ingoing and outgoing wave functions are computed as linear 
combinations of propagating plane waves at a given energy. The wave function in the GNR can 
be obtained using the quantum transmission boundary method 26 ' 27 . This allows us to calculate 
the spin- dependent transmission coefficients T±. Finally, we use the Landauer-Biittiker scattering 
formalism to calculate the spin-polarized electric current I± across the sample as follows 



where /(e) = [exp(e//cB^) + 1] 1 is the Fermi-Dirac distribution at temperature T (the value T = 4K 
was used thoughout the paper). 

SUPPLEMENT B. DIRAC THEORY FOR ARMCHAIR GNRS 

For energies close to the Dirac points, electrons in graphene can be described effectively by a 
two-dimensional two-valley Dirac equation , ! . The boundary conditions of GNRs require the wave 
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FIG. S4. a: Scheme of a GNR of width W = 8ao. The boundary conditions for wave functions can be obtained by adding two 
rows of atoms (plotted in gray) at y — 0, W + ao and setting the wavefunction to in those points, b: Transmission across a 
series of M potential steps. The incident plane wave with amplitude Ao splits into a reflected and a transmitted component 
with amplitudes Bo and Am+i, respectively. 

function to vanish on the (fictitious) sites just outside the GNR, i.e., at y = and y = W + ao, 
where the y axis is perpendicular to the direction of the GNR and the lower edge of the GNR is 
located at y = ao/2 (see Fig. S4 a). In the case of armchair GNRs, this affects both sublattices and 
the boundary conditions can be fulfilled by a superposition of two states from different valleys with 
the same energy E = hvF(k 2 _ + k 2 ) 1 ^ 2 and equal longitudinal momentum hk\\, but with opposite 
transverse momentum ±hk±, measured from the Dirac points 20 ' 28 . Here Vp is the Fermi velocity in 
graphene. Note that the effective description given by the Dirac equation holds as long as the k • p 
approximation remains valid, i.e., for not too narrow GNRs. 

Quantization of transverse momentum 

Since the valley momenta K and K' can be chosen parallel to k±, the transverse wave function 
can be written (p±(y) = sin [(If + k±)y\ where K = 47r/3ao- This function is evaluated on the 
honeycomb lattice with y e Nao/2 and oscillates rapidly. The transverse momentum k±, however, 
is small and quantized by the conditions (j)±(W + ao) = 0±(O) = 0. The allowed values for k± are 
given by (K + k± n )(W + a ) = nir, n G Z, and the spectrum reads 



E n {k\\) = ±hv F Jk 



h 2 



(S3) 



Taking into account that W is an integer multiple of a /2, one finds that the spectrum is gapless 
if 28 



W — (3ni + l)a /2, meN 



(S4) 



For asymmetric armchair GNR, as in Ref. 20, n\ is even. For symmetric armchair GNR, W is an 
integer multiple of ao and n\ is odd, such that W = (3n — l)a , n G N implies a gapless spectrum. 

In real samples there are small gaps even in the case (S4), which are due to edge effects 10,11 not 
included in the simple Dirac ansatz nor the homogeneous tight-binding formulation. In this work 
we consider symmetric armchair GNRs of width W = nao, where the integer n is different from 
3N — 1, e.g., W = 40a . These have a gap anyway and are quite robust against edge effects. Then, 
the allowed values of the transverse momentum are 



k 



3(W + a ) 

and the half-gap is E = Ei(0) = 7i-hv F /[3(W + a )]. 



n=l,2,4,5,7,8, 



(S5) 



Transfer-matrix description of transmission 



For geometries with potential depending only on the longitudinal coordinate x, the transverse 
momentum k± together with the wave function 4>±(y) is conserved, and it suffices to solve for the 
longitudinal wave function <p\\ (x). We consider the transmission across a piecewise constant potential 
profile, as sketched in Fig. S4 b. The solution of the Dirac equation for each spin o = ±1 and in 
each interval of constant potential value V is the superposition of two counter-propagating sublattice 
pseudo-spinors 
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with tan^ = ku/k± and kn = [(E — V) 2 / (hv-p) 2 — k"]] 1 ^ 2 . The solution may be evanescent because 
Eq. (S6) holds also for \E — V\ < hvp\k±\, when kn and 9 become imaginary. Then, the general form 
of the wave function in each slab j with potential Vj and momentum k\\ = kj is 

' ~ i " A ;{.r) A ,-(./•) 

(S7) 



e -i0j/2 e i0j/2 
_ e i0j/2 _ e -i6 j /2 



\ B A X ) ) 



y B A x ) ) 



where Aj(x) = Ae tkjX and Bj(x) = Be tkjX , such that 
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(S8) 



with L-j — Xj-\.\ — Xj. At each junction, kj changes but the wave function remains continuous 
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With the help of Eqs. (S8) and (S9), one writes down the transfer matrix for the whole system 



( A ^ 

A-M+l 



(S10) 



with d 



Q fl Q—l 



For the transmission problem depicted in Fig. S4 b, the boundary condition is no incoming electron 
from the right, Bm+i = 0. The reflection coefficient at the left is the ratio of reflected to incident 
current, R = \B \ 2 /\A \ 2 . For the transmission coefficient one has to take into account that the 
longitudinal momenta fcjw+i and ko are different if Vq ^ Vm+i, such that the ratio of transmitted to 
incident current is T = (\AM+i\ 2 kM+i)/(\Ao\ 2 ko). 

SUPPLEMENT C. BAND STRUCTURE OF THE UNBIASED LATTICE 



The Dirac formalism allows us to analytically study the system in the limit iV — > 00, when the 
energy regions with high transmission become transmissions bands surrounded by energy regions 
with T = 0. For an unbiased lattice with identical barriers of width a and spacing b, there are 
only two different transfer matrices involved, G a and Gt,. In the limit iV — > 00, the superlattice 
eigenf unctions have the Bloch phases exp(±iKl), that are the eigenvalues of the transfer matrix 
G — G a Gf) over one lattice period I — a + b. Thus, the dispersion relation E(K, k±) is obtained as 
cos(Kl) = Tr(G)/2, or again 22 ' 24 

1 



cos(Kl) = cos k a l a cos k b l b + 



. . sin k a l a sin k b l h . 
sin u n sin Vh 



(Sll) 



If I Tr((j)/2| > 1, there is no propagating solution with K real and E falls into the bandgap of the 

superlattice. In the main text, Fig. 2 a shows the bands of the infinite superlattice together with 

the transmission across a finite sample of five barriers. 
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